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Abstract. We revisit the problem of maximizing expected logarithmic utility 
from consumption over an infinite horizon in the Black-Scholes model with pro- 
portional transaction costs, as studied in the seminal paper of Davis and Nor- 
man [Math. Operation Research, 15, 1990]. Similarly to Kallsen and Muhlc- 
Karbe [Ann. Appl. Probab., 20, 2010], we tackle this problem by determining 
a shadow price, that is, a frictionless price process with values in the bid-ask 
spread which leads to the same optimization problem. However, we use a dif- 
ferent parametrization, which facilitates computation and verification. More- 
over, for small transaction costs, we determine fractional Taylor expansions of 
arbitrary order for the boundaries of the no-trade region and the value func- 
tion. This extends work of Janecek and Shreve [Finance Stoch., 8, 2004], who 
determined the leading terms of these power series. 



1. Introduction 

It is a classical problem of financial theory to maximize expected utility from 
consumption (cf., e.g., (TSJ [20] and the references therein). This is often called 
the Merton problem, because it was first formulated and solved in a continuous- 
time setting by Merton [TSJ QI5] . More specifically, he found that - for logarithmic 
or power utility and one risky asset following geometric Brownian motion - it is 
optimal to keep the fraction of wealth invested into stocks equal to a constant 9, 
which is known explicitly in terms of the model parameters. 

Magill and Constantinides [17] extended Merton's setting to incorporate propor- 
tional transaction costs. In particular, they showed that - again for logarithmic or 
power utility - it is optimal to engage in the minimal amount of trading necessary 
to keep the fraction of wealth in stocks inside some no-trade region [9, 9] around 9. 

Their somewhat heuristic derivation was made rigorous in the seminal paper 
of Davis and Norman [4] , who also showed how to compute 9, 9 by solving a free 
boundary problem. 
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Using the theory of viscosity solutions, this was further generalized by Shreve 
and Soner 23J. Moreover, the same approach allowed Janecek and Shreve [12] to 
derive rigorous first-order asymptotic expansions of 9, 9 and the value function for 
small transaction costs (also cf. (2[ [23j [21] [26] for related asymptotic results) . 

All of the above papers use methods from the theory of stochastic control. On 
the other hand, dating back to the pioneering papers of Jouini and Kallal [13] and 
Cvitanic and Karatzas [5] , much of the general theory for markets with transaction 
costs is formulated in terms of duality theory leading to consistent price systems 
resp. shadow prices. These are frictionless markets evolving within the bid- ask 
spread of the original market with transaction costs, which lead to an equivalent 
dual optimization problem. For logarithmic utility, Kallsen and Muhle-Karbe [T3] 
solved the Merton problem with transaction costs by simultaneously determining a 
shadow price and its optimal portfolio, i.e., by determining the dual and the primal 
optimizer at the same time (also compare |16) for an extension of this approach to 
a limit order market). 

In the present article, we revisit this problem using a different parametrization 
in the spirit of [B] . We employ arguments that are tailor-made for sufficiently small 
transaction costs, i.e., which focus on asymptotic expansions and thus seem promis- 
ing for tackling more complicated models. In particular, we determine rigorous 
asymptotic expansions both for the boundaries 9, 9 of the no-trade region and the 
corresponding value function, where terms of arbitrary order can be algorithmically 
computed. This extends the first-order expansions of j!2j . However, we emphasize 
that we only consider the particularly simple case of logarithmic utility, unlike [12] , 
who also deal with the more involved case of power utilities. An extension of the 
present approach to power utilities is subject of current research. 

The remainder of this article is organized as follows. After introducing the Mer- 
ton problem with transaction costs in Section 2, we present a heuristic derivation of 
our candidate shadow price process. Subsequently, in Section 4, we prove that this 
candidate is indeed well-defined and a shadow price process for sufficiently small 
transaction costs. Moreover, we derive asymptotic expansions for the boundaries of 
the no-trade region. Afterwards, we also determine an asymptotic expansion for the 
value function. Finally, in Section 6, we show that the shadow price corresponds 
to the solution of the dual problem. 

2. The Merton Problem with transaction costs 

We study the problem of maximizing expected logarithmic utility from consump- 
tion over an infinite horizon in the presence of proportional transaction costs, as in 
(HQ31E3]. Ffx a filtered probability space (Q, («^i)tem + ,P), whose filtration is 
generated by a standard Brownian motion (W t ) tg R + . We consider a market with 
two investment opportunities, namely a bond and a stock. The ask price process S 
of the stock is supposed to follow the Black-Scholes model 

(2.1) dSt/St = fjdt + adW t , S > 0, n,a > 0, 

and the bid price process is assumed to be given by (1 — X)S for some A 6 (0, 10 
This means that one has to pay the higher price St when purchasing the stock at 



This notation, also used in |25| . turns out to be convenient in the sequel. It is equivalent 
to the usual setup with the same constant proportional transaction costs for purchases and sales 
(compare, e.g., g][l2][l4][23]). Indeed, set S = and A = Then ((1 - X)S, S) coincides 
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time t, but only receives the lower price (1 — X)St when selling it. The price of the 
bond is supposed to be constant and equal to 5° = 1. As in |4| [14], we make the 
following standing assumption, which ensures that the holdings in bond and stock 
remain positive at all times. 

Standing Assumption 2.1. 

O<0:=4<1- 

The case 8 > 1 can be dealt with using minor modifications of the present 
approach (cf. [B]). The degenerate case 9 = 1 has to be treated separately, see 

[THUS]. 

Since transactions of infinite variation lead to immediate bankruptcy (compare 
PQ), we confine ourselves to the following set of trading strategies. 

Definition 2.2. A trading strategy is an M 2 -valued predictable process (tp°,tp) of 
finite variation. Here, (ipQ_,tpo-) = (rjB,rjs) G R^_\{0,0} represents the initial 
endowment in bonds and stocks, whereas tp® and tp t denote the units held in bond 
and in stock at time t, respectively. A (discounted) consumption rate is an R+- 
valued, adapted stochastic process k satisfying J K s ds < oo a.s. for all t > 0. A 
pair ((ip°,ip),K) of a trading strategy ((p°,ip) and a consumption rate n is called 
portfolio/consumption process. 

To capture the notion of a self-financing strategy, we use the intuition that no 
funds are added or withdrawn. As in [14], this leads to the following notion. 

Definition 2.3. A portfolio/consumption process ((tp° , tp), k) is called self-financing, 
if 

(2.2) tp° = tp° _ + f (1 - X)S t dtpj - f S t dipl - f K t dt, 

Jo Jo Jo 

where <p — tp' — ip^ for increasing predictable processes tp*, tp* which do not grow at 
the same time. 

Note that since S is continuous and tp is of finite variation, integration by parts 
yields that this definition coincides with the usual notion of self-financing strategies 
in the absence of transaction costs if we let A = 0. 

The subsequent definition requires the investor to be solvent at all times. For 
frictionless markets, i.e., if A = 0, this coincides with the usual notion of non- 
negativity of wealth. 

Definition 2.4. A self-financing portfolio / consumption process ((tp , ip), k) is called 
admissible, if the liquidation wealth process 

V t (<p°,cp) :=tp° t +tp+(l-X)St-tp t -S u t>0, 

of(tp ,tp) is a.s. nonnegative. 

We consider the following classical investment / consumption problem. 



with ((1 - X)S, (1 + X)S). Conversely, any bid-ask process ((1 - X)S, (1 + A)5) with A e (0, 1) 
equals ((1 - \)S, S) for S = (1 + X)S and A = -H.. 
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Definition 2.5 (Primal problem). Given an initial endowment {rjB,rjs), an ad- 
missible portfolio / consumption process ((ip , tp), k) is called optimal if it attains the 
maximum in 

(2.3) u (1 _ x)S!S (r)B,Vs) ■= sup E[f e~ st log( Kt )dt ) 

{{<p°,<p),k) \Jq J 

over all admissible portfolio/consumption processes ((<p , <p), k). Here S > denotes 
a fixed given impatience rate. 

Note that since 5 > 0, the value function of the Merton problem without trans- 
action costs is finite by [5J Theorem 2.1]. This upper bound for the value function 
in the present setup with transaction costs then implies that the latter is finite, 
since it is also bounded from below by the finite utilities that can be obtained from 
liquidating at time zero and consuming at a sufficiently small constant relative rate. 

The dual problem corresponding to (|2.3[) can be defined as follows (compare [3])- 

Definition 2.6 (Dual problem). Denote by U(t,y) = -e~ St (\og(e St y) + 1) the 
conjugate function of the utility function U(t,x) — e~~ St \og(x), i.e., the Legendre 
transform of —U(t, —x). Then a pair (S, Z) is called dual optimizer, if it attains 
the minimum in 



oc 



(2.4) «(i-a)s,s(»7b, Vs) := mf E I / U I t, ) dt ) + 

(S,z) \Jo V ${VB+VsS )J J 5 

Here (S, Z) runs through the set of all continuous semimartingales S with values 
in the bid-ask spread [(1 — A) 5, S] and the density processes Z of corresponding 
martingale measures, i.e., the positive martingales with initial value 1 for which 
SZ is a martingale. 

Remark 2.7. The pairs (S, Z) correspond to the consistent price systems of [91122], 
In the frictionless case A = there is only one of these, namely the stock price 
process S itself and the corresponding density process Z determined by Girsanov's 
theorem. In this case, the dual problem reduces to the usual notion for frictionless 
markets as in, e.g., [T51 Theorems 3.9.11 and 3.9.12]. The only difference is the 
constant 1/(5, which is added to let the optimal primal and dual values coincide 
(compare Theorem 16.11 below) . 

3. Heuristic derivation of the shadow price 

The key to solving both the primal and dual problems above is the following 
concept: 

Definition 3.1. A shadow price process S is a continuous semimartingale with 
values in the bid-ask spread [(1 — X)S, S], such that the maximal expected utilities 
U S s( 7 7B; 7 7S') * n th e frictionless market with price process S and un-\)s,s(VBiVs) 
in the original market with bid-ask process ((1 — A)^, S) coincide. 

From an economic point of view, it is obvious that the maximal expected utility 
for any frictionless price process S with values in the bid-ask spread [(1 — X)S, S] 
is at least as high as for the original problem with transaction costs, since one can 
always buy at least as cheaply and sell at least as expensively. To achieve equality 
of the maximal expected utilities, the optimal number tp t of stocks in the frictionless 
market with the shadow price process S should increase (resp. decrease) only when 
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St = St (resp. St = (1 — X)St). Otherwise one could achieve higher utility than in 
the original market with transaction costs. This portfolio will then also be feasible 
in the original market with transaction costs and hence be optimal for this market 
as well. 

It is well-known, dating back to |17) . that there is no trading in the original 
market with transaction costs as long as the fraction 

- - 1 

nt : " p° + <ptS t ~ 1 + tf/faSt) 

of total wealth in terms of the ask price S invested in stocks lies inside some 
interval [£, 6] around the Merton proportion 9 = n/cr 2 . Let us reparametrize this 
conveniently for the computations below (compare [6]): Set 

(3.1) c:=i^, s:=V-m 
1 ' 9 ' (1-6)0 

and define 

m t := ift/(c(pt). 

Then 7r G [9, 9] is equivalent to m £ [S/s, S] . We now make the ansatz that the 
shadow price S is given by 

S t = G(m t ,S t ), 

for some function G : R 2 — > R. This is equivalent to assuming that St is a function 
of the current ask price St and the current fraction ixt of wealth, in terms of the ask 
price S, invested in stocks. However, using mt instead of Ttt turns out to be more 
convenient in the calculations below. 

The function G enjoys a scaling property inherited from the logarithmic utility 
function and its constant relative risk aversion: 

G(amt,aSt) = aSt = aG(m t , St), Va > 0. 

This is because the optimal number ip® of bonds changes to atp® whereas the optimal 
number ip t of stocks and the dimensionless constant c remain the same when we 
rescale from St to aSt- Therefore, we can pass to the function 

g(S t /m t ) := G(l, S t /m t ), 

such that 

St = G(m t ,St) = m t g(St/mt). 

We now heuristically derive a free boundary problem that determines the func- 
tion g as well as the constants c and s. With these quantities at hand, we then 
explain how to construct the process m and in turn the shadow price S. 

If S is a shadow price, it is well-known that for logarithmic utility the optimal 
consumption rate at time t > in the frictionless market with price process S is 
given by 

K t = S(cp° + (p t St), 

where S > denotes the discount factor in the optimization problem (cf., e.g., [20)). 
Consequently, we have 

(3.2) dipt = 0, drf = -%t° + <pJh)dt = -6<p° t (l + g(St/m t )/c)dt 

in the no-trade region, where we have inserted our ansatz St = mtg(St/mt) and 
the definition mt = ipt/(apt)- Now consider the process 

Y t := St/mt = Stopt/ipt, 
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which represents the fraction of wealth in stocks (in terms of the ask price S) over 
wealth in bonds, normalized by c. In view of the dynamics (12. lj) of S and (|3.2|) of 
<p° and (f, Ito's formula yields 

(3.3) dYt/Yt = dS t /S t - d<p° t /tp° t = (fx+ 8(1 + g{Y t )/c))dt + adW t . 

Again by Ito's formula, the dynamics of the shadow price S = (S/Y)g{Y) are 

. 5/5 ( M'{Y t )Y t + £g"(Y t )Y? ~ S(g(Y t ) - g'(Y t )Y t )(l + g(Y t )/c) \ 
dS t /S t = ^ dt 



(og'{Y t )Y t 



dW t 



V 9(Yt) 
=: fi(Y t )dt + a(Y t )dW t , 

during each excursion into the no-trade region. Due to Morton's rule for logarithmic 
utility [T!5] , the Merton ratio for S must equal the proportion of total wealth invested 
in stocks in terms of the shadow price S, i.e., we must have 

HY t ) _ <ptSt g(Y) 

= 7T< := 



a(Y t )2 l - tf + nSt c + g(Y t y 
After rearranging, this yields the following ODE for g: 

(3.4) 9 » { y) = _ WM + + - ^M) . 

c + g{y) y °~ \ c j \ y z y ) 

Note that in the limit 6 — > 0, this is precisely the ODE from [B]. By ()3.1[) . the 
process Y = S/m takes values in [1,5]. Hence it is natural to complement the ODE 
()3.4|) with some boundary conditions at 1 and s as in [6j. Since we have a second 
order ODE, and the parameters c and s are also yet to be determined, we need two 
boundary conditions at each endpoint. The shadow price S — mg(S/m) should 
equal the bid price (1 — X)S at the selling boundary 9, where m = S/s, resp. the 
ask price S at the buying boundary 9, where m — S. Hence we impose 

(3.5) , 9 (1) = 1 and g(s) = (1 - A)5. 
Moreover, we add the smooth pasting conditions 

(3.6) </(!) = ! and g'(s) = 1 - A. 

They can be derived by looking at the ratio S/S, which should stay in [1 — A, 1], and 
setting its diffusion coefficient equal to zero at 1 — A and 1 (compare [5J). These four 
boundary conditions should now determine the function g as the solution to our 
second-order ODE as well as the two free variables c and s, which in turn identify 
the boundaries 9 and 9 of the no-trade region. 

With g,c,s at hand, it now remains to extend this construction from one excur- 
sion into the no-trade region to the entire positive real time line. To this end notice 
that, even when trading in the stock does happen, the process Y = aptSt/f® must 
always remain within [1, s], in order to keep the corresponding fraction of wealth in 
stocks in [9, 9}. However, the diffusion coefficient in the SDE (|3.3j) is bounded away 
from zero when Yt approaches either 1 or s. Therefore we need to complement the 
SDE (|3.3[) for Y with instantaneous reflection at both endpoints of [l,s]. Then 
Y will by definition remain in [l,s] and follow the claimed dynamics during each 
excursion into the no-trade region. With the process Y at hand, we can then define 
the candidate shadow price S = (S/Y)g(Y). 
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4. Construction of the shadow price 

We begin with an existence result for the free boundary value problem from the 
preceding section. As in [TH Proposition 4.2], it would be possible to obtain a 
similar result for any relative width A of the bid-ask spread. However, we restrict 
ourselves to the case of sufficiently small A here. This is all we need to derive asymp- 
totic expansions in terms of A below and seems much better suited to the analysis 
of more complicated models, since it requires only a much coarser inspection of the 
involved equations. 

Lemma 4.1. If A > is sufficiently small, then there are constants c> (l — 9)/6 
and s > 1 as well as a function g : [1, s] [1, (1 — X)s] satisfying the ODE 

( 4.D ,-(,) = -m» + *u+*M\m-m 

c+g(y) y a 2 \ c J \ y 2 y 
with boundary conditions 

(4.2) 9 (1) = 1, g{s) = {l-\)s, g'(l) = l, g'{s) = (1 - A). 

Proof. We first ignore the free boundary 5 and consider only the ODE (|4.1|) with 
initial conditions g(l) = g'(l) = 1, with c as free parameter. This problem can be 
written as 

d s g{s,c) = Fi(s,g,h,c), 
d s h{s,c) = F 2 (s,g,h,c), 
g{l,c) = h{l,c) = l, 

where F%(s, u, v, c) :— v, and 

„ , . 2v 2 20v 25 f u\ ( u v \ 

F 2 {s,u,v,c) :=— + — 1 + - — -- • 

c + u s a z \ c) \s z s/ 

Since c := {1 — 0)/ 9 { — 1, 0}, the function F 2 (and of course also Fx) is analytic in 
a complex neighbourhood of {s,u,v,c) — (1,1, l,c). By a standard existence and 
uniqueness result for analytic ODEs [1 1 } Theorem 1.1], we obtain a neighbourhood 
of (s,c) = (1, c) with a unique analytic solution g{s.c). The coefficients Oy of its 
Taylor expansion 



g{s, c) = 1 + (s - 1) + ^ J2 o« (a - lY{c - c? 

i>2 j>0 

can be calculated recursively. Note that a^j = a\j = Soj for j > 0, where S is the 
Kronecker delta, which follows from the initial conditions g{l,c) — g'{l,c) = 1. 

From this analytic function of two variables we will now construct the solution of 
the free boundary value problem (|4.1|) - (|4.2|) . by supplying appropriate quantities c 
and s (depending on A). These should satisfy 

(4.3) Xs = s-g{s, c) and g{s, c) - sd s \ s=s g{s, c) = 0. 

We divide the latter equation by s— 1, reflecting the fact that we are not interested in 
the solution s = 1 (which exists because of the initial conditions g{l) = g'{l) = 1). 
A Taylor expansion yields 

(4.4) SM-Sd^c) = E E bij(s 1)i(c _ 5)j 

i>Q j>0 
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for some real, explicitly computable coefficients bij. It turns out that 60,0 = 0, while 
^1,0 = 0(9 — 1) is non-zero, so that we can apply the implicit function theorem for 
multivariate analytic functions [101 Theorem LB. 4] to obtain an analytic function 

(4.5) H(c) = 1 + di(c - c) 1 = 1 + ^(c - c) + 0((c - cf) 

i>l 

such that (|4.4|) vanishes if H(c) is substituted for s. Here and in what follows, 
the coefficients of the power series symbolized by the 0(-)-s can be algorithmically 
computed. Now we insert the function H(c) into the right hand side of the first 
equation in (|4.3[) : 

H(c) - g(H(c),c) = a^(c - c) 3 + 0((c - c) 4 ), 

where do = 1(1 — 9) 2 /9 4 . If we insert H(c) also on the left hand side, the first 
equation in (|4.3[) becomes 

A(l + di(c - c) + 0((c - c) 2 )) = a^(c - c) 3 (l + 0(c - c)). 

Dividing by a^ 1 (1 + 0(c — c)) , we obtain 

\{a + O{c-c)) = (c-cf. 

Let us denote the coefficients of the series on the left hand side by a^: 

A(a + a x (c- c) + ...) = (c- c) 3 . 

Now take the third root: 

A 1 / 3 4 /3 (l + ^( c -c-) + ...)= c -c. 

(Note that the coefficients of powers of power series can be algorithmically calcu- 
lated [8].) Since ao = |(1 — 9) 2 /9 4 is non-zero, we can apply the implicit function 
theorem, so that the latter equation has a unique solution for small A. It is an 
analytic function of A 1 ' 3 : 

c-c = fl J /3 A 1 / 3 + 0(A 2 / 3 ). 

The coefficients in this power series expansion can be calculated with the Lagrange 
inversion formula: 

(4.6) [A fc / 3 ](c - c) = i[z k ^]a k /3 (l + + . . . )\ k>l, 

where the operator [z k ] extracts the fc-th coefficient of a power series. (This is 
analogous to the proof of Proposition 6.1 in [6].) With this value of c, the function 
g(s, c), together with s = H(c), is the desired solution of the free boundary value 
problem. To see that it maps [1, s] to [1, (1 — A)s], note that otherwise the equation 
g{s) — sg'(s) = would have a solution in the open interval ]1, s[. This is impossible, 
since H(c) is the unique zero of (|4.4[) . 

From (I4.6[) . we can calculate as many terms in the expansion of c as we want. 
The first of them are 
(4.7) 

l_± (^_\ 1/3 AV 3 + 2^ + 3(l- g )' / 6 \ ^ A2/3 + 



26» \9(l-9)J 129 \9(l-9) 
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Plugging this expansion into (|4.5|) gives an expansion of s, again to arbitrary order. 
The first terms are 

<«> - 1 + Gii^) 1/8 ai,s + K^>) 2 " +ow - 

Finally, note that (|4.7p shows that c > c = (1 — ff)/Q for sufficiently small A. □ 

Henceforth, we consider sufficiently small transaction costs A > for Lemma 
14. II to be valid and write g resp. c, s for the corresponding function resp. constants 
defined therein. 

The existence of the process Y as the solution to the SDE (|3.3p with reflection 
at 1 and s is a consequence of a classical result of Skorokhod [24] . 

Lemma 4.2. Let y G [l,s]. Then there exists a strong solution to the Skorokhod 
SDE 

(4.9) dY t =Y t (» + 5(l+g(Y t )/c))dt + Y t adW t , Y = y, 

with instantaneous reflection at 1 ands, that is, a continuous, adapted, [l,s]-valued 
process Y and nondecreasing adapted processes L and U increasing only on the sets 
{Y t = 1} and {Y t = s}, respectively, such that 

Y t = Y Q + f y„(/i + 6{1 + g{Y u )/c)du + f Y u adW u + f Y u dL u - f Y u dU u 



Jo Jo Jo Jo 

holds for all t £ K.+ . 

Remark 4.3. In the limit 5 — > 0, the process Y is just a geometric Brownian 
motion instantaneously reflected at 1 and s. Hence it coincides with the process 
S/m from |6]. Let us also point out that Y is related to the process (3 from [TU 
Lemma 4.3] by a nonlinear transformation. However, the present parametrization 
evidently leads to simpler dynamics. Moreover, the additional process C from [T4l 
Lemma 4.4] is no longer necessary. 



Proof of Lemma \4~J% Since g is bounded on [1,8] and c > (1 — 6)/ 6 > by 
Lemma 14.11 the coefficients of (|4.9I) are globally Lipschitz on [l,s]. Hence the 
claim follows from 24 . □ 

It was shown in [51 Proposition 2.1] that, for S — 0, the process S = mg(S/m) = 
(S/Y)g(Y) is an Ito process, i.e., does not involve local time. We now deduce 
from Lemma [4.21 that this still holds for arbitrary 5 > 0. The present proof only 
uses Ito's formula and the smooth pasting conditions (|4.2p for g. In particular, it 
allows to avoid the somewhat ad hoc localization procedures used in the proof of 
[5J Proposition 2.1]. 



Proposition 4.4. For y 6 [l,s] and Y as in Lemma \4-%\ S := (S/Y)g(Y) is a 
positive ltd process of the exponential form 

5 Sog(v) 



y 

with 



-£{X), 



dX t := dS t/ S t = I , dt + (M®gL\ dW t , X =0, 



(c + g(Y t ))g(Y t )J \ g(Y t ) 
and S takes values in the bid-ask spread [(1 — A)5, S 
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Proof. Using Ito's formula, we obtain 



dg(Y t ) = g'(Y t )Y t (fx + 5(1 + g(Y t )/c)dt + adW t + dL t - dU t ) + a — g" (Y t )Y? dt 



2 



d ^-.-o(l + ^)dt~dL t + dU t . 



and 

c 

" St/Yt 

Hence integration by parts yields 

,5,5 ( w'WYt + ^g"(Y t )Y? - S(g(Y t ) - g'(Y t )Y t )(l + g(Y t )/c) \ 
dS t /S t = ^ ^ j dt 

WW Y * dW t + (dL t - dU t ). 



V 9(Yt) J g(Y t ] 

Since L and U only increase on {Y t = 1} rcsp. {Y t — s}, where the integrand 
g'(Y t )Y t —g(Yt) is zero by the smooth pasting conditions for g, the last term vanishes. 
The Ito process decomposition of X now follows by inserting the ODE (|4.1[) for g" . 

For the second part of the assertion, notice that g"(l) < 0, because c > (l — ff)/6. 
Hence it follows from a comparison argument that g'(y)y~g{y) < on [1, s]. As this 
implies that the derivative (g'(y)y — g(y))/y 2 of g(y)/y is non-positive, <?(1)/1 = 1 
and g(s) /s = 1 - A yield that S = Sg(Y) JY is indeed [(1 - X)S, S]-valued. □ 

Using standard results for frictionless markets, we can now determine the optimal 
portfolio/consumption pair for S = (S/Y)g(Y). 

Lemma 4.5. Set 

!1, if TIB > Socqs, 

s, ifrj B <S cris/s, 

Socrjs/f]B otherwise. 



For Y and S as in Lemma \4.2\ and Proposition ^. 4\ respectively, define 

1 



as well as 



l + c/g(Y t )' 

V t = {vb + Vs§ )£ ( [ ^dS t - [ Sdt 
\Jo b t Jo 

Kt = SVt, 

<Pt=*tV t /St, ip° t =(l-n t )V t . 



Then 



(a ii\ , f* VuCg'(Y u )Y u 1T [ f y u cg>(Y u )Y u 

Jo 9(Yu){c + g(Y u )) J g(Y u )(c + g(Y u )) 

and ((tp , <f), k) is an optimal portfolio/consumption pair for the frictionless market 
with price process S . The corresponding wealth process and the optimal fraction of 
wealth invested in stocks in terms of S are given by V and fr, respectively. 
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Note that the first resp. second cases in ()4.10|) occur if the initial fraction of 
wealth in stock r]sSo/(r]B+VsSo) lies below 8 — (1+c) -1 resp. above 8 = (l+c/s)" 1 . 
In this case, there is a jump from the initial position (</?[]_, V?0-) = (vb, Vs)i which 
moves the initial fraction of wealth to the nearest boundary of the interval [8_, 8] . 
Since this initial bulk trade involves the puchase resp. sale of stocks, the initial 
value of S is chosen so as to match the initial ask resp. bid price in this case. 

Proof. Since dS/S — dX t :— fi{Y t )dt + a(Y t )dW t is an ltd process with bounded 
coefficients and a is bounded away from zero, its optimal portfolio/consumption pair 
is characterized by standard results for frictionless markets, cf., e.g., [3 Theorem 
3.1]. In particular, the optimal fraction of wealth (in terms of S) invested in stocks 
is given by Merton's rule, i.e., equals 

fi{Y t ) _ g(Y t ) _ . 
a(Y t ) 2 c + g(Y t ) ** 

as claimed. Moreover, the optimal wealth process, the optimal consumption rate 
and the optimal numbers of bonds resp. stocks are indeed given by V, k, (p a , and (p. 
Now Ito's formula and the ODE (|4.1[) for g imply 

\c + g{Y t )J \{c + g(Yt)yj \{c + g(Yt)) 2 

Moreover, by Yor's formula, we have 
£(J £^dX t -J Q 5dt) 

= S {\ , C fv , d((X,X) t -X t )- f Sdt 
\Jo c + g(Y t ) J 

, , c 2 a 2 g'(Y t ) 2 Y 2 f c f cag'(Y t )Y t 



/o g(Y t ) 2 (c + g(Y t )) 2 J Q J Q g(Y t )(c + g(Y t )) 

Taking into account that the process L — U is continuous and of finite variation, 
integration by parts finally yields the claimed representation for tp. □ 

The representation (|4.11j) and the definition of So via (|4.10|) imply that stocks 
are only purchased (resp. sold) when S = S (resp. S = (1 — X)S). Hence the 
optimal portfolio/consumption process ((<p ,<p),K) for S is also admissible for the 
bid-ask process ((1 — X)S,S). Since shares can be bought at least as cheaply and 
sold as least as expensively when trading S instead of (1 — X)S,S), it follows that 
the maximal expected utilities coincide, i.e., S is a shadow price. Made precise, 
this is the content of the following analogue of [131 Theorem 4.6]. 

Theorem 4.6. The portfolio/consumption process [{<p ,<p),k) from Lemma \4.5\ is 
also optimal for the bid- ask process ((1 — X)S, S). In particular, S is a shadow price 
in this market. 

Proof. This follows verbatim as in the proof of [TU Theorem 4.6]. □ 

In view of Lemma l4 . 5 1 and the definition of g, the optimal proportion tt of stocks in 
terms of the shadow price S takes values in the interval [(l+c) -1 , (l+c/((l— A)s) -1 ]. 
Translating this into terms of the ask price S, we find that the optimal fraction 
7r in terms of S is kept within [(1 + c) _1 , (1 + c/s) -1 ], i.e., the lower resp. upper 
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boundaries of the no-trade region are given by 6 — (1 + c)~ 4 resp. = (1 + c/s) -1 . 
The following series expansions are immediate consequences of (|4.7p and (|4.8|) (upon 
taking one additional term in both expansions). 

Corollary 4.7. The lower and upper boundaries of the no-trade region in terms of 
the ask price S have the expansions 

9 = J—=9-(-0 2 (l- 9)A A 1 / 3 -^L( —^—] A 2 / 3 

0(9(1 - 0) 2 (1 - 20 + 20 2 )- 4<5 2 /cr 4 + 65/^(3 - 80 + 50 2 )) 6 

—7 — A 



360(1 - 0) 



+ 0(A 4 / 3 ) 

and 

1 



0(3(1 - 0) 2 (3 - 266» + 260 2 ) - ^/a 4 + 6<5/cr 2 (3 - 80 + 50 2 )) 6 



360(1 -0) 0(1 -0) 

+ 0(A 4 / 3 ), 

respectively. Consequently, the size of the no-trade region in terms of S satisfies 



l + c/s 1 + c 
: (60 2 (1 - 0) 2 ) 1 / 3 A 1 / 3 

0(3(1 - 0) 2 (3 - 160 + 160 2 ) - 4£ 2 Ar 4 + 65/a 2 (3 - 89 + 50 2 )) 6 

180(1-0) 0(1-0) 

0(A 4 / 3 ). 



A 



Further terms of these expansions can be algorithmically computed, if desired. 

The first correction terms are of order A 1 / 3 , which has been conjectured in [23j 
Remark B.3] and proved in [12]. Whereas the first-order corrections are symmetric 
around the Merton proportion 0, the second correction terms of order A 2 / 3 are both 
negative. This means that it is optimal to start transacting for smaller fractions 
of stock than indicated by the first-order expansions. As pointed out by |12) . 
an intuitive explanation is that consumption from the bank account continuously 
increases the proportion of stock. Accordingly, the A 2 / 3 -terms vanish in the limit 
for 6 — > 0. Interestingly, however, the impatience rate 8 > only shows up starting 
from the third term of the expansion for the width of the no-trade region even in 
the presence of consumption. 



Remark 4.8. We now compare the expansions from Corollary 14.71 to the results 
of Janecek and Shreve |12j mentioned above. We rename their parameter A as A; 
it is related to our A by A = A/ (2 — A) (see the footnote in Section^. The bounds 
of the no-trade region become 
(4.12) 

1 /3 

1 =0- f-0 2 (l-0) 2 ) A 1 / 3 -^(120 2 (l-0) 2 )- 1 /3p/3 + o( ^ ) 



l + 2c/(2-A) V 2 
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and 
(4.13) 

1 „ /3„n,, „o\ ;,/, 298 



002(1-0)2^ A 1 /3_^(1202( 1 _ )2 ) -l/3^2/3 +O( ^ 



1 + 2c/((2 - X)s) 

The A 1 / 3 -terms agree with those established in [TJ], while our A 2 / 3 -terms are not 
those that were conjectured in |12j . based on heuristic considerations. 

We are indebted to Steve Shreve, who kindly checked these calculations again. 
While the arguments worked out in [TJ] turned out to be perfectly correct, he found 
that an assumption on which these calculations were based did not hold true. It 
was stated in equation (4.27) of 12 that "there is considerable evidence" that the 
coefficient of the A term in the expansion of the value function has a certain nonzero 
value, while subsequent calculations by Steve Shreve indicate that this coefficient 
is zero. When this term is set to zero, then heuristic calculations based on the 
approach of [H] lead to (|4.12l) and (|4.13p as well. 

5. Characterization and asymptotic expansion of the value function 

In the frictionless market with price process S, standard control methods allow 
to characterize the value function for the primal problem from Definition 12.41 as 
the unique classical solution to the corresponding HJB equation. This yields the 
maximal utility in the frictionless shadow market, which coincides with the maximal 
utility in the original market with bid-ask process ((1 — X)S, S) by Theorem[ 



Theorem 5.1. Let y, Y, and S be defined as in Lemma \4-5\ Then, if A > is 
sufficiently small, we have 

(5-1) u {1 ^x)S,s{vb,vs) = log(?7B + VsSo)/S + w(y), 

where w is the unique (classical) solution to the ODE 

(5.2) ?-s 2 w"(s) + ( fi + S ( 1 + ^] ) sw'(s) - 8w(s) + jfj?) . + log(<5) -1=0, 



2 w V V c / 7 28d 2 {s) 
on [l,s] with Neumann boundary conditions 

w'{\) = w'(s) = 0. 



Proof. We first show that the boundary value problem (|5.2p has a unique solution, 
which is an analytic function. This is similar to the first part of the proof of 
Lemma I4TT1 We first solve the ODE (|5.2I) with initial conditions w(l) = Wi+£ and 
w'(l) = 0, with arbitrary £, and where 

H 2 -25a 2 log 6 

is the value to which w(l) converges for A — > 0. This yields an analytic function 
w(s,c,0 = w 1 +^ + J2 E a l]k {s-\ye{c-c)\ 

i>2 j,k>0 

for certain real coefficients a^. (Recall that c tends to c = (1 — 9)/ 8 as A — » 0.) 
To satisfy the boundary condition w'(s) = 0, we solve the equation 



= 

s—s(c) 
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for £. After inserting the expansion ()4.5p for s, we obtain a solution £(c), analytic 
at c = c. Then u>(s) = w(s, £(c(A)), c(A)) solves the boundary value problem. 

Now fix A > sufficiently small, such that Lemma [4.1[ Corollary 14.71 and the 
above existence result for w can be applied. Then, in particular, the mapping 
y h- > can be extended to a C 2 -mapping on the whole real line, e.g., by attaching 
suitable parabolas at y = 1 and y = s. Ito's formula in turn shows that 



M t := f e- Su \og( Ku )du + e' St 
Jo 



/o 
--M 



log(Vt) 



w(Y t ) 



t e-^(^^+aY u w'(Y u ))dW u , 



where we have inserted the dynamics of V (cf. Lemma 4.5) and Y (cf. Lemma 
4.2) as well as the ODE and boundary conditions for w. In particular, the process 
(M t )t>o is a local martingale. By the proof of Proposition 14. 4\ we have a(y) = 
<rg'{y)/(9(y)/y) < <V(1 - A) for y e [1, s], such that < ff < (1 + c/((l - A)s)) _1 < 
(1 + c/s) _1 and 1 < y < s imply 



E 



-2Sv 



d-(Y u )jr v 



aY u w'(Y u \ 



1_ / cr/(l - A) 
" 2(5 U(l + c/a) 



(TStO 



for |«/|max : = max j/G[i,s] l w '(2/)l- Therefore the process (M{)t>o is in fact a true 
martingale and also uniformly integrable, since it is bounded in L 2 (P) by the Ito 
isometry. Hence M t converges in i 1 (P) to a random variable M M with 

(5.3) E{M 00 ) = B(M ) = log(Vb)/<5 + w(Y ) = log(ry s + vs^/S + w(y). 

Now notice that e~ St w(Y t ) — > 0, because w(-) is bounded on [l,s]. Moreover, as 
log(Vt) is an Ito process with bounded coefficients, it follows that e~ St log(Vt)/<5 — > 
in L 2 (P) and hence also in probability. Therefore Mt — ► Jq 00 e~ St log(Kt)dt in 
probability, such that the a.s. uniqueness of limits in probability yields 



E 



,-<5t 1 



log(/st)dt = J B(M 00 ). 



Combined with f|5 . 3[) . this completes the proof. 



□ 



It does not seem possible to come up with an explicit solution to the above ODE 
even in terms of the function g, unlike in the case without consumption treated in 
[51 125) . However, we can again derive an expansion in powers of A 1 / 3 . Note that 
this expansion cannot be recovered by letting the relative risk aversion tend to 1 in 
the corresponding result for power utility [12] . The reason is that the value function 
for power utility does not converge to the value function for log utility. 

Corollary 5.2. In addition to the assumptions of Theorem I5.il suppose that our 
portfolio starts on the Merton line, i.e., 

VsSo _ ft 
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Then the value function satisfies 

U(i-\)s,s{vb,Vs) =^ log(5(?7B + vsSq)) + ^ (j^ - 5 

(5 4) 3^( M (^-^)V3 

(5 ' 4) _ 4-21/3^10/3 +C; W' 

T/iis expansion can be continued to any order, with algorithmically computable 
coefficients. 

Proof. We expand the first term in (|5.ip : 

- log(77 B + 7756*0) = g log(?7s + VsS Q g{y)/y) 

= i logfe + r; s 5 ) + ± log + 8 - 1 

= ^log(T7 B + r) S S ) + 0(A). 

The last line follows after expanding log into a Taylor series, and inserting y = c/c 
and the series for g and c from the proof of Lemma 14.11 As for the second term 
in ()5.1|) , we use the series for w found in the proof of Theorem 15.11 and obtain 



() 1 /V aU 10 ^ 3 2 / 3 ( M (a 2 - M )) 4 / 3 x2/3 

= V 2^ - * J + 4-2i/3^io/3 + ° (A) - 

Adding both expansions yields (|5.4p . □ 

Note that the expansion (|5.4p can easily be refined to a bivariate expansion in 
terms of A 1 / 3 and 7r — 9, valid for A > sufficiently small and n G [0, 6*] inside the 
no-trade region. Its first terms arc identical to (|5.4[) ; ttq — 8 only appears from the 
0(A 4 / 3 )-term onwards. 

6. Solution of the dual problem 

We conclude by showing that the shadow price S and the density process of 
the corresponding martingale measure are a solution of the dual problem from 
Definition [ 



Theorem 6.1. Let Y and S be defined as Lemma \4.5[ and set 

/o °Vt) 
Then 

U(l-\)S,s(llB,Vs) = V(i-\)s,s(VB, Vs), 
i.e., there is no duality gap in the market with bid-ask process ((1 — \)S,S), and 
(S, Z) is a dual optimizer. 

Proof. First notice that for any continuous semimartingale S with values in the 
bid-ask spread [(1 — X)S, S], it follows as in the proof P3J Theorem 4.6] that 

U(l-X)S,s(VB,Vs) < us,s(vb,Vs)- 
Together with, e.g., [7J Lemma 2.3], this implies 

U(1-X)S,S(VB,VS) < V(i_x) StS (VB,Vs)- 
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For the shadow price S, we have vn_\\s,s(VB, Vs) < v § s( 7 7S; 7 7s) by definition and 
U S s^BtVs) = u (i-\)s,s{' r lB,Vs) by Theorem 14.61 It therefore remains to show 

v S,s(Vb,Vs) = u ss (t]b,T]s), 

i.e., that there is no duality gap in the frictionless market with price process S. As 
the frictionless market with price process S is standard and complete in the sense 
of [15, Definition 1.7.3], this follows from }15^ Theorem 3.9.12] and we are done. □ 
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